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Abstract 

The scope of constrained differential renormalization is to provide renor¬ 
malized expressions for Feynman graphs, preserving at the same time the 
Ward identities of the theory. It has been shown recently that this can be 
done consistently at least to one loop for abelian and non-abelian gauge theo¬ 
ries. We briefly review these results, evaluate as an example the gluon selfen¬ 
ergy in both coordinate and momentum space, and comment on anomalies. 


1 Introduction 

Differential renormalization (DR) is a renormalization method in coordinate 
space which provides finite Green functions without any intermediate reg¬ 
ulator or counterterms Q. It has been applied to a variety of problems (a 
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complete list of references can be found in Ref. [0). As initially formulated, 
DR is too general for it introduces many arbitrary constants in the process 
of renormalization, which have to be fixed at the end of the calculations 
requiring the fulfilment of the relevant Ward identities. A constrained ver¬ 
sion has been recently proposed which aims to fix this arbitrariness from the 
beginning and automatically fulfil the Ward identities ill- In Section 2 
we briefly review this constrained procedure. It has been explicitly shown 
to work at one loop for gauge theories (QED, Scalar QED and QCD are 
discussed in Refs. i i, i and ||, respectively). The method is defined 
in coordinate space, but one can perform the calculations also in momen¬ 
tum space, using the Fourier transforms of the renormalized expressions in 
coordinate space. A computer code doing all operations automatically is 
available ||. For illustration, we evaluate in Section 3 the one-loop gluon 
selfenergy in both coordinate and momentum space, and verify that it is 
transverse, as required by gauge invariance 0. Finally we comment on 
anomalies in Section 4 and on future prospects in the conclusions. 

2 Overview of constrained differential renormal¬ 
ization 

Constrained differential renormalization (CDR) is essentially a set of rules 
which allow to expand any Feynman graph in coordinate space in a set of ba¬ 
sic functions and at the same time determine the (universal) renormalization 
of these functions. The renormalized graph is then defined as the correspond¬ 
ing linear combination of the renormalized basic functions. There are two 
types of rules, those stating how to manipulate the singular diagrams and 
those fixing the renormalization of the basic functions. The former establish 
that all the algebra (including Dirac algebra in four dimensions) must be 
simplified first, treating the singular expressions as if they were well-defined. 
This includes contracting all possible Lorentz indices and using the Leibniz 
rule for derivatives. After this, the Feynman graph is written as a linear 
combination of total (external) derivatives of basic functions. These are 
products of propagators with all (internal) derivatives acting only on the 
last propagator: 

A mi (zi)A m2 (z 2 ) ■ ■ ■ 0^ zl A mn (zi + z 2 -\ -1- z n - 1 ), (1) 
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A m [l] 

A m [c^] 

Bmp 777-2 [1] 

B777,!7772 [ d / l ] 

B777,!777,2 P] 
Bt77!777,2 [<9^5^] 


T777,1777,2777-3 P] 

^777,-1777,2777.3 

Tttt,!777,27773 [^<9/7 ] 

Tttt!777,27773 



Q?77-i 7772m37774 

Q?77l777277737774 P^/7^] 
Qt77 1777277737774 [ d /^ Oi / dpdcr ] 





Table 1: Singular basic functions for renormalizable gauge theories in four 
dimensions in the Feynman gauge. Lines are ordered according to the num¬ 
ber of propagators and columns according to the degree of singularity. The 
function A, that appears in tadpoles, is defined as A m [0] = OA m (z)5(z). 

where A mi (zi) = _L^ ; w ith K\ a modified Bessel function ||, is 

the propagator of a particle of mass ] and 0^ Zl is a differential operator 
of order r acting on the variable z\. Zi = Xi — Xi + 1 , i = 1 ,.. . , n — 1 are 
the coordinate differences and we work in four-dimensional euclidean space. 
The rules of the second type reduce the degree of singularity (changing sin¬ 
gular expressions by derivatives of less singular ones, which are prescribed to 
act formally by parts on test functions) and extend mathematical identities 
among tempered distributions to arbitrary functions. In four-dimensional 
renormalizable theories (in the Feynman gauge) there are only thirteen sin¬ 
gular basic functions. They can have one (A), two (B), three (T) or four 
(Q) propagators, corresponding to the lines in Table 1, and the singular be¬ 
haviour at coincident points can be of order z~ 4 (logarithmic), (linear), 
z~ 6 (quadratic) and z~ 7 (cubic), corresponding to the columns in that ta¬ 
ble. A cubic singular behaviour only occurs in tadpole diagrams that vanish 
for symmetry reasons in usual theories. A detailed description of the CDR 
rules and the renormalized expressions of these singular basic functions in 
both coordinate and momentum space are given in Ref. ||]. With those 
expressions and the decomposition described above one can systematically 
calculate any one-loop amplitude. A computer program performing all op¬ 
erations automatically is available. Its description can be found in Ref ||]. 

1 The expression of the propagator reduces to A (zt) = 4 ^ 3-77 in the massless case. 
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3 The gluon selfenergy 


To show how the method works, let us evaluate the gluon selfenergy in 
QCD. We shall first do the calculation in coordinate space, and afterwards 
perform a Fourier transform to obtain the corresponding momentum space 
expression. To simplify the resulting expressions, we shall consider the case 
of massless quarks, although the massive case is analogous. The lagrangian 
and Feynman rules in coordinate space can be found in Ref. The four 



Figure 1: Feynman diagrams contributing to the gluon selfenergy. 


contributing diagrams are depicted in Fig. 1. Their “bare” expressions in 
coordinate space, in the Feynman gauge, are 


yrab ( A ) 


TTtlfe (B) 

TTClb (C) 

TTtlf) (D) 


l -g 2 N c 5 ab A(x) 

bp,o{Dp — dp) + 5op(dp— dp,) + b P p(do ~D a ) 

fiav(dp — Dp) + b V p{D a — do) + bpo{dv —d u ) 

A{x) , (2) 

-6g 2 N c 5 ab 5p U A(x)5(x ), 

(3) 

-g 2 N c 5 ab dpA(x)d l/ A(x ), 

(4) 

-\9 2 N f 5 ab Tr[^A(x)7 M ^A(x)7, / ] , 

(5) 


where x = z\ = x\ — X 2 is the coordinate difference, N c the number of colours 
and Nj the number of quark flavors. The partial derivatives d a = act 
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on the next propagator, while d a act on the previous one. D a are partial 
derivatives acting on the external gluons. Using the Leibniz rule, they reduce 
to minus total derivatives of the amputated expression. We first expand 
Eqs. (ID in basic functions: 


n“t (j4) = \g 2 N c d ab {(2d ll d u -55^n)B[l]+5(d fl B[d u ]+d u B[d „]) 

+ 2ViB[9 ff ] - 10B[d,A] - 2<5 #1I/ B[m]} , (6) 

n“t (B) = -6 g 2 N c 5 ab S^A[l], (7) 

K {C) = -9 2 N c 5 ab [d,B[d u ] - Bld^])} , (8) 

n“t (D) = -2 g 2 N f 5 ab {d^[d u \ + d v B[d»} - S^B^} - 2B[d fl d u ] 

+ <5/ii/B[m]} . (9) 


Then, we replace the basic functions by their renormalized expressions in 


Table 2 of Ref. 

namely 


A r [1] = 

0, 

(10) 

B r [1] = 

1 ^logx 2 -/!/ 2 

(ID 

647T 4 X 2 

B R [^] = 

\d,B R [l], 

(12) 

B R [a} = 

o, 

(13) 

B^d,] = 

3 (d/A ^ I1U D)B R [1] + 2887 ^ 2 {d^d u 5^U)5{x) 

,(14) 


where the derivatives are prescribed to act formally by parts and M is the 
renormalization scale. The result for the gluon selfenergy is the sum of 
Eqs. (§ D with A[l] and B \0\ replaced by A R [1] and B n [O ]: 


it 


ab R 
fiv 


i 


1447T 2 


g 2 6 ab (d„d v - S^a) 


1 lo£x 2 -/\L 2 

(15W - 6N f )^n^- 2 -+ (2 N c - 2 N f )S(x) 


(15) 


The result is transverse = 0), as required by gauge invariance @. 

is less singular than the bare 11“^ and admits a finite Fourier trans¬ 
form in four dimensions: 

K R (P) = J d 4 xe“-m;t fl (x). (16) 


5 









This momentum space gluon selfenergy is easily obtained acting with the 
external derivatives by parts on the exponential, as prescribed by DR, and 
Fourier transforming log ^7 [|J and <5(x): 


■prafe R 


1447T 2 


g 2 S ab {p^Pu - S^p 2 ) 
M 2 


M 

(15 N c - 6 N f ) log — r - 2N C + 2N f 


(17) 


where M = 2M/^e and 7 e = 1.781... is Euler’s constant. In more involved 
cases including masses one can use the Fourier transforms of basic functions 
in Appendix B of Ref. Q. 

Alternatively, one can do all the computations in momentum space, once 
the Fourier transforms of basic functions are known. Every step in the 
calculation in coordinate space can be translated into momentum space. 


4 Anomalies 

When dealing with singular expressions, anomalies can appear at the quan¬ 
tum level. In some cases (spurious anomalies) the symmetry may be re¬ 
stored with finite local counterterms, while in others the anomaly cannot 
be avoided and has physical relevance. This is the case of the ABJ trian¬ 
gular anomaly ||. In general, the source of anomalies in CDR is the non¬ 
commutation of renormalization with the contraction of Lorentz indices, 
which comes from the fact that the rules we have imposed are incompatible 
with this operation a 11 - Let us briefly discuss the case of the ABJ 
anomaly in QED. In the one-loop triangular diagram of one axial and two 
vector currents, the trace of six Dirac matrices and one 75 appears. Al¬ 
though this trace is unique in four dimensions, it can be written in different 
ways that lead to different expansions in (singular) basic functions (each of 
them with different contractions of the internal derivatives). The final re¬ 
sults for different expansions differ by finite pieces, and in all cases at least 
one of the three independent Ward identities is broken. We have checked 
that the democratic form, 

Tr[757A7a7z/767/i7c] = 4 {e Xa vb^c - av^bc + ^Xavc^ 

“ 1 “ €\ab/j,3 vc 4 “ ^Xvbc^afi ^Xabc^/-LV 4 “ ^Xafic^vb ^Xvfic^ba 4 “ bfic^av 

4 “ ^avbfi^Xc 4 “ ^av/ic^Xb ^avbc^fiX 4 “ ^-vbfic^Xa 4 “ ^-vb/iX^ca 
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preserves the vector Ward identities, giving the correct anomaly in the axial 
current [|]. A detailed discussion of anomalies within constrained differential 
renormalization will be presented elsewhere. 


5 Conclusions 

Constrained differential renormalization is a new version of DR which ren¬ 
ders renormalized amplitudes fulfilling automatically the Ward identities 
of symmetric theories. It has been explicitly shown to work at one loop 
for renormalizable gauge theories. Here we have verified that the one-loop 
gluon selfenergy renormalized with this method satisfies the transversality 
condition implied by gauge invariance, and have briefly discussed how the 
ABJ triangular anomaly is recovered. One must also wonder about higher 
order calculations. They require generalizing the renormalization rules (in 
particular, the treatment of integrated internal points must be specified) as 
well as the set of basic functions. This issue is under study. 
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